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ABSTRACT 



Context. Using the technique of spectral disentangling, it is possible to determine the individual spectra of the components of a mul- 
tiple star system from composite spectra observed at a range of orbital phases. This method has several advantages: it is unaffected by 
line blending, does not use template spectra, and returns individual component spectra with very high signal-to-noise ratios. 
Aims. The disentangled spectra of a binary star system are very well suited to spectroscopic analysis but for one problem: the absolute 
spectral line depths are unknown because this information is not contained in the original spectra (unless there is one taken in eclipse) 
without making assumptions about the spectral characteristics of the component stars. Here we present a method for obtaining the 
atmospheric parameters of the component stars by the constrained fitting of synthetic spectra to observed and disentangled spectra. 
Methods. Disentangled spectra are fitted using synthetic spectra and a genetic algorithm in order to determine the effective tempera- 
tures, surface gravities and relative light contributions of the two stars in a binary system. The method is demonstrated on synthetic 
spectra and then applied to the eclipsing binary V615 Per, a member of the young open cluster NGC 869 (hPersei). 
Results. The method works well for disentangled spectra with signal-to-noise ratios of 100 or more. For V615 Per we find a normal 
He abundance but an Mg abundance, which indicates bulk metallicity, a factor of two lower than typical for nearby OB stars. 
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1. Introduction 

The technique of spectral disentangling (spd) allows the isola- 
tion of the individual spectra of the component stars of a double- 
lined spectroscopic binary system from a set of composite spec- 
tra observed over a range of orbital phas es. It w as original ly for- 
mulated in the wavelength domain b y Sim on & Sturm d 19941) 
and in the Fourier domain bv lHadraval (Il995h . The techiiique si- 
multaneously returns the best-fitting individual spectra and the 
orbital velocity amplit udes of the two stars. A detaile d overview 
of SPD can be found in lPavlovski & Hensbergd (l2009l) . 

Compared to other methods of radial velocity measurement, 
' SPD has several advantages. Firstly, it is independent of template 
' spectra so avoids any systematic eiTors due to spectral differ- 
ences between the target and template stars. Secondly, it is not 
' affected by the b lendin g of spectral lines of th e two stars (see 
! lHensbergeet"al] 120001: ISouthworth & Claused l2007h . Thirdly, 
the resulting disenta ngled spectra contain the comb ined signal 
of the input spectra ('Pavl ovski & South worthll2009h so have a 
much higher signal-to-noise (S/N) ratio. 

There are two disadvantages of the spd approach. The first 
of these is that the continuum normalisation of the input spectra 
has to be very good in order to avoid low-fre quency spurious 
patte rns in the resulting disentangled spectra dHensberge et alJ 
I200&) . The second is that relative continuum light contributions 
of the two stars cannot be found using spd as this information 
is itself not contained in the observed spectra, unless a spectrum 
has been obtained during an eclipse (Iliiic et al. 2004). 

SPD is well suited to the spectral analysis of stars in binary 
systems. Each disentangled spectrum contains only features due 



to one star, so can be analysed using standard methods for single 
stars. The high S/N ratios of disentangled spectra are very help- 
ful to this process, but the undetermined continuum light ratio 
between the component stars complicates the spectral analysis. 
In this work we present a method to fit synthetic spectra to dis- 
entangled spectra, where the atmospheric parameters of the stars 
are determined simultaneously with the relative light contribu- 
tions of the stars. A genetic algorithm is used for the optimiza- 
tion in order to ensure that the best solution is found in a param- 
eter space which suffers from strong degeneracies, in particular 
between effective temperature (Teff) and surface gravity (logg). 

An important application of spd is the study of detached 
eclipsing binary star systems (dEBs). These represent the pri- 
mary source of directly-measured masses an d radii of stars , 
and as such are cornerstones of stellar physics ( Andersenl l 19911 : 
ToiTes et al. 2010). spd can be used to measure the velocity am- 
plitudes of the stars, which are necessary for the mass and ra- 
dius measurements, simultaneously with obtaining the individ- 
ual stellar spectra for spectral ana lysis (Pavlovski & Hensbergd 
20051: iPavlovski et al. '2009^. iPavlovski & South worth 2009). A 
major advantage of dEBs to this process is that the surface 
gravities of the stars can be obtained to within +0.01 dex from 
the mass and radius measurements: these parameters can then 
be fixed in the spectral an alysis and thus the degeneracy be- 
tween T^g and log g avoided dSimon et al.ll994l : lHensberge et al.l 
2000). 

In this work we demonstrate the genetic algorithm approach 
to fitting disentangled spectra on the dEB V615Persei, a mem- 
ber of the young open cluster h Persei. ISouthworth et alJ (l2004aL 
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Fig. 1. Synthesized Hy line profiles (red lines) and the best- 
fitting synthetic spectra (blue lines) for the case with S/N =100 
in Table[T] The primary star is below the secondary star on the 
plot. 



hereafter SMS04) obtained a series of high-resolution spec- 
tra of V615 Per a nd analysed them with published light curves 
dKrzesinski et al.l[l999.) to measure the masses (4.08 and 3.18 
Mo) to accuracies of 2% and the radii (2.29 and 1.90 Rq) to 5%, 
resulting in surface gravities measured to within 0.05 dex. The 
Teff values were found to be 15 000 + 500 K and 11 000 ± 500 K. 
SMS04 found that stellar evolutionary models required a sub- 
solar metal abundance (Z 0.01) to reproduce the measured 
masses and radii of V615 Per. 

The Perseus Double Cluster comprises h Persei (NGC 869) 
and;\f Persei (NGC 884). It has been extensively studied via deep 
CCD photome t ry ( Keller et al. 2001; Marco & Bernabeu 2001; 
Slesnick et all l2002i: iCapilla & Fabregat ,2002: .Currie et aL 



2O10l) . from which there is general agreement on its distance 
(2.3 to 2.4 kpc) and age (13-14Myr). But these studies assumed 
a solar chemical composition, and their results may be system- 
atically wrong if this assumption is incorrect. 

Conflicting results on the chemical composition of the 
Perseus Double Cluster are present in the literature. Detailed 
abund ance analyses based on high-resolution spectra of hot 
stars dLennon et alJ 119881: iDufton et alJ [1990) hav e chal- 
lenge d previous findings of l ow he l ium abundances jNissenl 
1976 1 'Klochkova & Panchuk" '1987'; 'Wolff^ & Heaslevl Il985h . 
Dufton et al. ( 1990) and Smartt & Rolleston ( 1997) found an ap- 
proximately solar metal abundance from high-re solution spec- 
tra, but this was not supported bv Vrancken et al.l 12000). In this 
work we attempt to shed additional light on this subject by mea- 
suring the helium and metal abundances of the stars in the dEB 
V615Per. 



2. Constrained fitting of disentangled spectra using 
a genetic algorithm 

A computer code has been constructed which fits synthetic spec- 
tra to disentangled spectra of a binary system in order to de- 
termine the atmospheric parameters Teif , log g, projected rota- 
tional velocities v sin /, Doppler shifts, and the light factors. The 
light factors are an important part of the analysis, and are pa- 
rameterised as LF, the fraction of the total system light pro- 
duced by one star for the wavelength or wavelength interval un- 
der consideration. The LFs for the binary components should 
sum to unity. The determination of the atmospheric parame- 



ters represents a difficult optimization problem, for which we 
use a genetic algorithm to minimise the of the fit to the 
data (Holland 1975; Charbonneau 1995). Our implementation 
is called genfitt (oENet i c FiiT ing) and in approach is similar 
to that of iMokiem et al.l (l2005h . Error estimates come from the 
covariance matrix, which is constructed using the Levenberg- 
Marquardt method. 

In order to save computing time we pre-calculate grids of 
synthetic spectra, genfitt linearly interpolates between these in 
Teft and log g, and then convolves them with a rotational profile 
using the rotin.f code of I. HubenjQ. The LTE grid covers T^s 
from 6000 to 15 000 K and the non-LTE grid covers Teff from 
14000 to 35 000 K. Both grids contain \ogg values of 2.5 to 5.0 
(cgs), and are stepped by 250 K in Tetf and 0. 1 dex in log g. 

There can be strong degeneracies between the fitted atmo- 
spheric parameters, most notably Teff and log^ for Balmer line 
profiles. This degeneracy between Teff and \ogg can be avoided 
by analysing dEBs, because their surface gravities can be known 
to within 0.01 dex from measurements of their masses and radii. 
In many cases the light ratio of the stars in a dEB can be obtained 
from the light curve analysis (e.g. Southworth et al. 2004b), so 
the LFs can then be fixed to known valuefl 

In order to test the performance of genfitt we synthesized 
disentangled spectra covering the HyS, Hy and H<5 lines, using 
representative atmospheric parameters and with Gaussian noise 
added to produce S/N ratios ranging from 25 to infinity. We then 
used genfitt to fit all the atmospheric parameters to the synthe- 
sized spectra, with the only constraint that the two light factors 
sum to unity. The results are given in Table[T] and show that for 
high S/N ratios (>100) the genfitt results reproduce the input 
Teif and log g values satisfactorily. For lower S/N ratios the in- 
herent degeneracy of these parameters causes them both to be 
underestimated by our method. By contrast, the LFs are repro- 
duced to well within the errorbars for all S/N ratios considered. 
An example fit is shown in Fig.[T] 

We obtained a second set of solutions in an 'unconstrained 
mode' where the LFs were not required to sum to unity. The 
results were, as expected, similar to but slightly poorer than the 
'constrained mode'. Finally, a third set of solutions were made 
with fixed log g values (Tabled, as would often be the case when 
analysing a dEB. We find that the situation is similar to that for 
the first set of solutions. The main limitation on the quality of 
these results is the degeneracy between T^s and log g: we find 
a correlation coefficient of 0.98 between these parameters for 
both components. The correlation with the LFs is much weaker, 
which is why the LF values are reliable even for low-S/N spectra. 

3. Application to V61 5 Persei 

genfitt has already be en used for studying the dEBs V380Cyg 
dPavlovski et alj|200 9') and V621 Per (Southworth et al. 201 1, in 
prep.). In both of these cases spectra of a high S/N ratio were 
available and genfitt returned excellent results. Here we chal- 
lenge it with spectra of a much lower S/N ratio. 

25 spectra of the dEB V615Per were obtained by SMS04, 
with a reciprocal dispersion of 0.11 ApixeL', a resolution of 
0.2 A and an average S/N of w50. They cover 4220^500A so 
include Hy, a number of helium lines, and the Mg II 448 1 A dou- 
blet and are well distributed through one orbital cycle, spd was 



' http : //nova . astro . umd . edu/index . html 

^ The light ratio can be poorly determined in some dEBs, and in this 
case obtaining a spectroscopic light ratio is an important part of m od- 
elling the light curves. For an example see lSouthworth et d1 l l2007h . 
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Table 1. Results of fitting synthesized disentangled spectra with genfitt for different S/N values. The output parameters and 
uncertainties are given, along with the difference compared to the input parameters and the reduced (xl)- 



S/N 


Teffl (K) 


ATeff, 


loggi 


Aloggi 


LFi 


A LFi 


T,ff2 (K) 


ATeff. 


logg2 


Alogg2 


LF2 


ALF2 




oo 


17991 


-9 


3.495 


-0.005 


0.775 


-0.005 


11 005 


-5 


3.995 


-0.005 


0.225 


+0.005 






±56 




±0.016 




±0.012 




±56 




±0.016 




±0.012 






500 


17 896 


-104 


3.492 


-0.008 


0.771 


-0.009 


10 896 


-103 


3.992 


-0.008 


0.229 


+0.009 


1.054 




+ 153 




±0.025 




±0.024 




±157 




±0.023 




±0.021 






200 


17 835 


-165 


3.497 


-0.003 


0.771 


-0.009 


10837 


-164 


3.987 


-0.013 


0.229 


+0.009 


1.032 




+ 198 




±0.027 




±0.015 




±211 




±0.024 




±0.019 






100 


17 753 


-247 


3.485 


-0.015 


0.769 


-0.011 


10748 


-252 


3.986 


-0.015 


0.231 


+0.011 


0.841 




±365 




±0.033 




±0.023 




±358 




±0.036 




±0.021 






50 


17 041 


-959 


3.421 


-0.079 


0.761 


-0.019 


10030 


-970 


3.929 


-0.071 


0.239 


+0.019 


0.774 




±441 




±0.035 




±0.081 




±430 




±0.036 




±0.086 







Notes. The input atmospheric parameters are: Teffi = 18000 K, loggj = 3.5, vsini'i = 105 kms LFi = 0.78, Tj-n^ = 11 000 K, loggj = 4.0, 
V sin (2 = 95 km s"' , LF2 = 0.22 

Table 2. Same as Table[T]but for fits with \ogg values for the two stars fixed. 



S/N 


Teffl (K) 


ATeffl 


LF, 


ALFi 


TeffJ (K) 


Areff2 


LF2 


ALF2 


xi 


00 


17 921 


-79 


0.776 


+0.001 


10921 


-79 


0.224 


-0.001 






±109 




±0.010 




±115 




±0.014 






500 


17 897 


-103 


0.771 


-0.004 


10909 


-91 


0.229 


+0.009 


1.065 




±160 




±0.022 




±166 




±0.021 






200 


17 825 


-175 


0.769 


-0.011 


10 808 


-192 


0.231 


+0.011 


1.042 




±236 




±0.018 




±236 




±0.018 






100 


17 750 


-250 


0.767 


-0.013 


10738 


-262 


0.233 


+0.013 


0.848 




±352 




±0.027 




±343 




±0.025 






50 


17 056 


-944 


0.756 


-0.024 


10030 


-970 


0.244 


+0.024 


0.777 




±416 




±0.082 




±422 




±0.085 







perfo rmed in Fourier space using the fdbinar^ code (llliiic et al.l 
I2OO4I) . The disentangled spectra have S/N values of about 160 
for the primary star (star A) and 80 for the secondary (starB). 

Since the spectra of V615Per cover only a limited wave- 
length range, our estimate of the TeffS of the component stars 
is restricted to the Hy line. Helium lines are also good Teft in- 
dicators, but instead we will use these later to obtain the helium 
abundance of the binary. Because the available light curves of 
V615Per are not definitive, the logg values for the two stars 
are known to modest accuracies of 0.059 and 0.050 dex. We 
therefore included them as genfitt fitted parameters (solution 
A), along with T^ff. We also obtained a solution B for compari- 
son, where the log g values were fixed to those found by SMS04. 
The LFs were constrained to sum to unity. For v sin / we adopted 
28 + 5 and 8 + 5 kms"' plus an instrumental broadening of 
16 kms ' (SMS04). The results are given in Tabled 

For the chemical abundance analysis we adopted the atmo- 
spheric parameters from solution A. Syn thetic spectra were cal- 
culated using ATLAs9 model atmospheres (Kurucz'lQTQ*) and non- 
LTE theoretical line profiles from the detail and surface codes 
(iGiddingsll 1 98 it iButler 1984j). A canonical microtu rbulence ve- 
locity of 2 kms"' was adopted (iTrundle et al.ll2007l) as we have 
too few spectral lines to fit for it.. The helium abundance of 
star A was derived via minimisation between the observed 
profiles of He I 4388 A and 4471 A, and profiles calculated for 
abundances in the range e(He) = 0.06-0.15. The results for both 
lines are given in Tableland correspond to a mean helium abun- 
dance of e(He) - 0.091 + 007. This abundance is solar to within 



^ |http : // sail . fer . zep/fdbinary/| 



Table 4. Chemical abundances derived for the components of 
V615Per. 



Component 


Species 


Wavelength (A) 


Abundance 


Star A 


He I 


4388 


0.090 ± 0.005 


Star A 


He I 


4471 


0.092 ± 0.005 


Star A 


He I 


mean 


0.091 ±0.007 


Star A 


Mg II 


4481 


7.26 ± 0.03 


StarB 


Mg II 


4481 


7.16 ±0.06 



the uncertainty [eo(He) — 0.089], so no deviations in helium 
abundance are detected for V6 15 Per A. The helium lines from 
StarB are too weak to be useful (Fig.|2]l. 

Magnesium is an excellent metallicity indicator for B-type 
stars as it does not particip ate in the CN O process so is un- 
modified by stellar evolution jLvubimkov e t al. 2005). These au- 
thors found a mean abundance of loge(Mg) = 7.59 + 0.15 for 
nearby B stars with reliable microturbulence velocities, in excel- 
lent agreement with the solar value of log eo(Mg) = 7.55 ± 0.02. 
Mg 11 448 1 A is a prominent feature in the spectra of both com- 
ponents of V 15 Per Theoretical line profiles were calculated in 
non-LTE for star A and in LTE for starB and the abundances ob- 
tained hy minimisation (Fig.|2|i. We find a mean abundance of 
log e(Mg) - 7. 21 + 0.07. This is 0.46 dex lower than the mean 
value found by Lvubimkov et al. (2005) and 0. 16 dex lower than 
that found by Daflon et al. (2003) for OB stars in the solar circle. 
The components of V615 Per have a subsolar Mg abu ndance and 
resemble halo B stars more than nearby examples dPaflon et al.l 
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Table 3. Results from genfitt analysis of the Hy line for V615 Per, with a comparison to the values found by SMS04. 



Parameter 


This work (solution A) 
Star A Star B 


This work (solution B) 
Star A Star B 


SMS04 
Star A Star B 


Teff (K) 

logg [cgs] 
LF 


14 710 + 210 11 520 ±290 
4.302 ± 0.035 4.361 ± 0.030 
0.676 ± 0.006 0.324 ± 0.006 


14920+ 190 11420 + 250 
4.328 fixed 4.381 fixed 
0.677 ± 0.004 0.323 ± 0.004 


15 000 + 500 11000 + 500 
4.328 + 0.059 4.381 +0.050 
0.65 + 0.03 0.35 + 0.03 




4460 4465 4470 4475 4480 



Fig. 2. Comparison between disentangled and renormalised 
spectra of star A (upper panel) and star B (bottom panel) with 
synthetic spectra calculated for the He and Mg abundances given 
in Table|4](blue lines). Only the Mgll lines and the Hei line for 
star A are fitted here. Other lines are shown but not used in our 
results. 

l2004l) . This implies that the hPersei open cluster has a subsolar 
metal abundance. 



4. Summary 

Spectral disentangling is a method for obtaining the individual 
spectra of the components of a binary star system from compos- 
ite spectra obtained at a range of orbital phases. A disadvantage 
of this method is that the continuum light ratios of the stars are 
not found, because this information is not present in the observed 
spectra without making assumptions about the spectral charac- 
teristics of the stars. We present the genfitt program, which uses 
a genetic algorithm to fit synthetic spectra to the disentangled 
spectra of both components of a binary system simultaneously. 
It returns the best-fitting atmospheric parameters (refr and log g) 
and the light contributions of the two stars. From tests with syn- 
thesized spectra we find that genfitt performs extremely well in 
determining the light ratio. It also returns reliable Teff and log g 
values in those cases where S /N of the input disentangled spec- 
tra is >100, which is the usual situation for observational studies. 

The light contributions of the two stars will normally 
sum to unity, which provides a useful constraint for genfitt. 
Contaminating light from a third star can in principle be found, 
in cases when the light contributions of the two stars in the bi- 
nary sum to less than unity. Once the light contributions of the 
stars have been found, their disentangled spectra can be renor- 
malised to the correct continuum levels. The resulting spectra 
can then be analysed using standard methods for single stars. 
If the two stars are eclipsing, their surface gravity values may 



be found to high precision and accuracy from analysis of the 
orbital velocity amplitudes found by spectral disentangling and 
light curves covering the eclipses. 

As a demonstration of the method we applied genfitt to 
spectra of the eclipsing system V615Per, a member of the 
hPersei open cluster. The metal abundance of this cluster is 
controversial (see Sect.[T]i but important in measuring its dis- 
tance by the isochrone method. The spectra were disentangled 
and fed into genfitt, and an abundance analysis was performed 
on the resulting renormalised spectra. The atmospheric param- 
eters returned by genfitt are in good agreement with previous 
work (SMS04) but are more precise. We find a normal solar he- 
lium abundance for V6 15 Per A (starB is cooler and has only 
weak helium lines). The magnesium abundances for both stars 
are lower than those found for nearby OB stars, indicating that 
h Persei has a subsolar metallicity. This is in agreement with the 
results of SMS04, based on the complimentary method of com- 
paring the masses and radii of the stars to the predictions of the- 
oretical stellar evolutionary models. 
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